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We show how an adequate post-Newtonian generahzation can be obtained for Newtonian dark 
matter halos associated with an empiric density profile. Applying this approach to halos that follow 
from the well known numerical simulations of Navarro, Frenk and White (NFW), we derive all 
dynamical variables and show that NFW halos approximately follow an ideal gas type of equation 
of state which fits very well to a polytropic relation in the region outside the core. This fact suggests 
that "outer" regions of NFW halos might be related to equilibrium states in the non-extensive 
Statistical Mechanics formalism proposed by Tsallis. 



PACS numbers: 04.20.-q, 02.40.-k 



I. INTRODUCTION 



The issue of dark matter clustering in halos of virial- 
ized galactic structures is one of the most interestingopen 
problems in astrophysics and cosmology [H 0, |3j M 1^ ■ 
The physical properties of this dark matter are uncer- 
tain, leading to various proposed physical matter models: 
thermal sources, meaning a colissionless gas of weakly in- 
teracting massive particles (WIMP's), which can be very 
massive (m ~ 100 — 200 GeV) supersymmetric ("cold 
dark matter" CDM) or self-interacting less massive [m ~ 
KeV) particles dQ ( "warm" DM). Otherproposals in- 
clude scalar fields (real and complex) |^, global mo- 
mopoles T^, axions, etc. However, all these models must 
comply from inferred direct and indirect observations 
that reveal the presence of DM: velocity profiles of rotat- 
ing stars, microlensing and tidal effects affecting satelite 
galaxies and galaxies within galactic clusters. Galactic 
DM is mixed with visible baryonic matter (stars and gas) 
clustering in galactic disks, making up about 5-10 % of 
the total galactic mass. Hence, it is a good approxima- 
tion to identify the gravitational field of a galaxy with 
that of its DM halo, considering visible matter as "test 
particles" in this field [ll|l3- 

Realistic galactic halos are obviously not spherically 
symmetric, but they are approximatelly so, since their 
global rotation is not dynamically significant Hence, 
we will consider throughout this paper that halos are 
spherically symmetric equilibrium configurations. As- 
suming the CDM paradigm and spherical symmetry, we 
can distinguish two types of halo models: those obtained 
from a Kinetic Theory approach, whether based on spe- 
cific theoretical considerations or on convenient ansatzes 
that fix a distribution function satisfying Vlassov's equa- 
tion [T^ , or those emerging from "universal" density pro- 
files obtained empirically by N-body numerical simmu- 
lations 15, 16, ITJ- In this paper we will consider the 
latter apprach, based on the well known numerical sim- 
ulations of Navarro, Frenk and White (NFW) ITsf. 



Although these simulations yield virialized equilibrium 
structures that reasonably fit CDM structures at a cos- 
mological scale (> 100 Mpc), some of their predictions in 
smaller scales ("cuspy" density profiles and excess sub- 
structure) seem to be at odds with observations 19, "2^, 
especially those based on galaxies with low surface bright- 
ness (LSB), which are supposed to be overwhelmingly 
dominated by DM and so well suited to examine the pre- 
dictions of various DM models. 

Galactic halos are newtonian systems characterized by 
typical velocities, ranging from 5 — 10 km/sec for dwarf 
galaxies up to about 1000 — 3000 km/sec for rich clus- 
ters. However, we believe that a study of these systems 
under General Relativity, as a post-Newtonian approxi- 
mation, might provide new information that can be use- 
ful and interesting for gravitational lensing and for the 
study of the interplay between cosmological scale evolu- 
tion and galactic DM. In any case, since General Rela- 
tivity is the best available theory of classical gravity, it 
is relevant from a theoretical point of view to be able 
to construct spacetimes that are suitable for important 
self-gravitating structures like galactic halos. 

All DM halo models derive a full set of dynamical vari- 
ables from a given "mass-density" profile. In a post- 
newtonian relativistic generalization we will assume that 
this density is the dominant rest-mass contribution to 
the matter-energy density, made up by rest-mass and 
an internal energy term proportional to suitably de- 
fined temperature and pressure (of a kinetic nature). 
Thus, we will assume an "ideal gas" type of equation 
of state [2^ in which this internal energy density 

becomes determined only by the hydrostatic equations 
themselves in the case of isotropic velocity distributions. 
In the anistropic case, which we leave for a future pa- 
per a, various empirical ansatzes can be assumed in or- 
der to relate radial and angular components of the stress 
tensor. 
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II. NEWTONIAN NFW GALACTIC HALOS. 

Following the CDM paradigm and assuming spherical 
symmetry, galactic halos must satisfy the following New- 
tonian equations of hydrostatic equilibrium 



M ' = 4 TT p 

, AttGM 
= ^ — , 



as well as the Navier-Stokes equation 

, 2a 
P' = P, 

r 

where P is the "radial" pressure and 

P-P\ 



P 



(1) 
(2) 

(3) 

(4) 



is the anisotropy factor relating P with the tangential 
pressure P^. Since we have three equations for five un- 
knowns [p, M, $, P, a), this system can be made deter- 
mined if two of these five functions becomes specified, for 
example, by assuming an "equation of state" somehow 
relating P and Pi with p. In the case of N-body numer- 
ical simmulations, we have virialized structures whose 
density profile p = p{r) can b e app roximately fit to a 
"universal" empirical function |l3l lisl [l8l | . The New- 
tonian system becomes determined once we have 
this density profile together with a suitable expression for 
a. In general, the simmulations yield anisotropic velocity 
distributions, so that specific ansatzes can be assumed or 
prescribed [23 for a 0. 

The well known N-body numerical simmulations by 
Navarro, Frenk and White (NFW) yield the following 
"universal" expression for the density profile of virialized 
galactic halo structures 0, : 



where 



5n = 



X {l + x) 



Ac3 



2 ' 



Po 



3[ln(l + co) -co/(l + co)]' 
Mo 



Pcrit = 253.8 /i^r^o 



kpc"^ 



X ^ —, 
r.. 



Co 



(5) 

(6) 

(7) 

(8) 



The virial radius r^i, is given in terms of the virial mass 
Mvir by the condition that average halo density equals A 
times the cosmological density po .2] 



4 7rM, 
3r3 



(9) 



where A is a model-dependent numerical factor (for a 
ACDM model with total f^o = 1 we have A ~ 100 



The concentration 
terms of M^,^ by [2 



arameter Cq can be expressed in 



Co 



62.1 X 



(10) 



where —0.5 < e < 0.5. Hence all quantities depend on a 
single free parameter A/^j, with a dispersion range given 
by e for different halo concentrations. The NFW mass 
function and Newtonian potential follow from integrating 
(HJ and I© for p given by © 



47rrf So Po 



ln(l + x)- 



1 + x 



ln(l + x) 



complying with the boundary conditions 

M«p„(0) = 0, M«p„(r_) = Af_, 
-$NFw(0) = = 47rG(5o/9orf, 



(11) 
(12) 



(13) 
(14) 



while circular rotation velocity (normalized by the char- 
acteristic velocity Vq) is simply 



47rG'A/ 



^0 



ln(l -t-x) 



1 



1 + a; 



(15) 



Given Pnfw, Af^pw and $nfw, radial and tangential pres- 
sures follow from integrating Q for a specific choice of 
a. There are analytic solutions of © for a = (isotropic 
case) and for various empiric expressions for a. A tho- 
rugh Newtonian treatement of NFW halos is found in 

Notice that according to the density profile (|3J) we have 
a diverging density at the symmetry center (x = 0). This 
is obviously an unphysical feature and points out to the 
fact that (jSJ has not been derived from any theoretical 
argumentation, but is simply a convenient empirical for- 
mula that fits the outcome of the NFW numerical sim- 
ulations which show that p^pw cx l/x near the central 
halo region. Since the virial radius r^j^ is the physical 
size of the resulting halos and these simulations are un- 
able to provide adequate resolution for distances to the 
halo center smaller than approximately 1% of the virial 
radius 01 ; all halo quantities presented here are, strictly 
speaking, only valid between a minimal r ~ 0.01 r^ir and 
Tvir (see section VII for further discussion on this issue). 



III. THE SPACETIMES OF GALACTIC HALOS 

As a good approximation p9| , we can consider the met- 
ric of a galactic halo to be given by a suitable "weak field" 
limit of the spherically symmetric static line element 



dr^ 



1 - KoM/r' 



% (16) 
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where <J>(r) is the relativistic generahzation of the New- 
tonian gravitational potential and kq — 2G/c^, so that 
M(r) has mass units. We will assume a momentum- 
energy tensor of the form 



T'^''^ fiu'^u'' +ph'''' 



(17) 



where u'^ = exp(-$/c2) h'''' = g''^ + u" and H"'' 
is the anisotropic and traceless (H" a = 0) stress tensor, 
which for the metric (|16|l takes the form 



n" b = diag [0, -2n, n, n] 



(18) 



so that p and 11 relate to P and Pi by 



Fl - p = 3 n, 



2Pi 



P = 3p. 



The field equations and momentum balance (T"'' .f, 
associated with ((TC|) - ((T^ are 



M' 
P' 



kqC^ M + 47rPr3/c2 
^ ~2 r (r - Ko Af) ' 

= -(^ + P) — P, 



(19) 
0) 

(20) 
(21) 

(22) 



where a is given by I0J. These equations are the rel- 
ativistic generalization of Q-®, though we must pro- 
vide a relation between /i and p. Since the particles in 
the collisionless gas making up galactic halos are interact- 
ing very weakly and the velocity anisotropics tend to be 
small: < a < 0.2 it is reasonable to assume that 
it is nearly an ideal gas and that total matter-energy 
density, /i, is the sum of a dominant contribution from 
rest-mass density, p (? , and an internal energy term that 
is proportional to the pressure P and to the velocity dis- 
persion = l^v^) ~ (w^). Hence it is reasonable to 
assume the equation of state of a non-relativistic (but 
non-Newtonian) ideal gas |2H 



p = pc 



. 3 a2 



P = pa\ 



(23) 



where the velocity dispersion is related to a kinetic tem- 
perature T by 



^2 P KT 

a = — ^ , 

p m 



(24) 



Since characteristic velocities in galactic halos are New- 
tonian, we have a^jf? ^ 1 and p k, pt? and so P ~ 
Pj^ ^ pc^, so that (|23|) provides a plausible equation of 
state for a relativistic generalization of galactic halos. It 
is evident also that in the newtonian limit j (? ^ 1 we 
recover the equilibrium equations Q-jSJ. 

What needs to be done now is to insert the equation 
of state (|^ into the field equations ipni) - (|^ . which 
becomes a set of equations that, just like be- 
comes determined once we specify the functional relation 
p = p(r) or p = p(*^') ^iid Q;(r). However, we do not 



need all the three equations H20() - (|22|l . since numerical 
simulations yield the density profile Q, we will assume 



p(t) 



Pnfw 



(25) 



and eliminate $' from (|21(l and l|22(l . leading to the fol- 
lowing set: 



M' = 4 7r 



3 P 
2 ? 



(26) 



P' = - 



r (r — Kq M) 



r 



(27) 



which becomes fully determined once we know a{r). We 
will solve these equations in a post-Newtonian approx- 
imation by keeping only terms up to order jc?. The 
velocity dispersion a (and/or T) can be obtained after- 
wards from P through and if^ . 

As mentioned before, we have 0. A 

quick calculation of the Ricci scalar using H21|l . I|26(l and 

m 



R 



SttG 



Pnfw C 



-I 



p 



(28) 



reveals the existence of a curvature singularity as r 0. 
This is wholy unphysical, since NFW halos within a gen- 
eral relativistic treatment should be weak field static 
spacetimes. However, as mentioned in the last parra- 
graph of section II, all variables associated with NFW 
simulations are valid in the range 0.01 < r/r^^^ < 1. In 
section VII we discuss how this issue can be dealt with 
appropriately. 



IV. POST-NEWTONIAN EQUATIONS FOR 
NFW HALOS. 



It is convenient to work with the following adimen- 
sional variables 



Y 
M 
V 



1 



M 



r.3 



eg AM 



AttSq p^rl 3 (5o M^^, '' 
P 



5q po 



(29) 
(30) 
(31) 



where the structural parameters Sq, po, co, A and Vq 
have been introduced in section II. The field equations 
(I21|-l|ll|) now becomes 



dM 

da; 

dP 

da; 



Y + ^eV x\ (32) 

[y+feP] [M + ePx^^] 2a ^ 
X [a; - 2 e M] x 

(33) 
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where 



ANALYTIC SOLUTIONS 



^0 



(34) 



so that in the hmit e ^ we recover the Newtonian 
equations ((TJ, ^ and The system (|22J|-lI23I) can be 
integrated by demanding that A4 and V comply with ap- 
propriate boundary and initial conditions. Since we have 
to use the explicit form of Y in (|29|l . then the analytic 
or numerical solutions of H32|l - H33() for specific choices 
of a, boundary conditions depend on M^j, through the 
definitions ® and ifTUI) . 

The metric function M = Vq rg M follows from (|32|l , 
while $ can be obtained by integrating 



_d_ 

dy 



M + eVy^ 
y[y-2eMy 



(35) 



The relativistic generalization of the Newtonian rotation 
velocity profile are the velocities of test observers along 
circular geodesies. From 0, 0, these velocities are 
= r which in terms of the adimensional variables 
becomes 



M+eVy^ 
y-2eM 



(36) 



Since Vq for typical galactic halos ranges from a few 
km/sec to ^ 3000 km/sec, the post-Newtonian cor- 
rections of order Vq/c^ will be very small: between 
0(e) ~ IQ-^ and 0{e) ~ 10"^. The post-Newtonian 
system associated with (|32|l - (|33() can be given as 



dM 
dx 

dP 

dx 

while 1)35(1 and H36|l become 



Yx^ + 0{e), 
Y M 2a 



V + Oie), 



d 

dy 



vi 



^ + 0(e), 
x 



(37) 
(38) 

(39) 
(40) 



In the post-Newtonian equations given above, V is de- 
cupled from A4 and <&, thus (irrespective of the form of 
V) the metric elements for all NFW halo spacetimes are 
up to order e 



- gtt 

Qrr 



,^^_M^^^Oie% (41) 



1 



2eM 



ln(l 



1 



- + 0(e2),(42) 

X \ + x\ 

where the metric functions M and 4> obtained from H37() 
and compl y w ith the boundary conditions ((T^ and 
lO (see also Notice that M and $ are finite 

at the center, even if Y diverges. Also, even if M di- 
verges as r — > oo, the metric components shown above 
are well behaved asymptotically, tending to flat space- 
time: —gtt 1 and g^r —'■1 as a; — ^ oo. However, 
the NFW spacetimes do not comply with the regularity 
at the center of a spherically symetric spacetime which 
requires the vanishing of all spacelike gradients (such as 
M' and $'). In fact, the Ricci scalar (|28|l in the post- 
Newtonain limit becomes: R = —2Ye + 0{e'^), hence 
there is a curvature singularity in the center even if the 
metric functions do not diverge. 



Even if all NFW halos have the same rest-mass density 
Y and metric functions A/, $, the form for the pressure 
depends on the assumptions one might make about a and 
suitable boundary conditions. For the remaining of this 
paper we will consider only the case of isotropic velocity 
distributions, leading to a = 0. In this case, the post- 
Newtonian Navier-Stokes equation ()38|l has the following 
analytic solution: 



P ^ [ln(l + x)f + A{x) ln(l + a;) - ^ In a; + 3 dilog (1 + x) - B{x) + 0(e), 



A{x) 



l-ix 



bx"^ 



2x2(1 



B{x) 



1 + 9a; 4- Ix^ 
2x(l + x)2 



(43) 



where C is a constant and the dilogarithmic function is In order to determine C, we need to examine the bound- 
defined as ary conditions of V. 

Intdt 



dilog(?/) 



5 



VI. POLYTROPIC EQUATION OF STATE 

The asymptotic behavior (x ^ 1) of P in (|43|l 

impHes that an asymptoticaUy flat configuration arises if 
we choose C ~ 7r^/2, so that V 0, scahng asymptot- 
icaUy as P oc In x/x"^. Since Y scales asymptoticafly as 
1/x^, this indicates a sort of power law relation between 
V and Y that (at least asymptotically) might be similar 
to a polytropic relation of the form 

r^KY^+^/", (45) 

where K and n (polytropic index) are constants. In order 
to examine the functional relation between Y and V, we 
provide in figure ^ a logarithmic plot ofV vsY (or equiv- 
alently InP vs InpVjj^), for the asymptotically flat case 
with C ^ applied to a halo with M,i, = 10^^ Mg, 

corresponding to a virial radius marked by a; = cq ^ 15 
(~ 150 kpc). For theoretical reference we show the curve 
associated with a polytropic relation (|45|) with n « 5.5 
and K « exp(— 1.7). As shown by the figure, the asymp- 
totically fiat NFW configuration fits very well this poly- 
trope, except for high density values corresponding to 
smaller x, up to the value x = xq ~ 0.01 cq that marks 
the resolution limit of numerical simulations (^ 1 kpc). 
This behavior is reasonable, since closer to the center 
{x close to xo) the NFW density profile becomes cuspy, 
while polytropic density profiles are characterized by a 
"fiat core" [ij. 



InV 
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FIG. 1: The thick curve is the plot of in "P vs in F, equiva- 
lent to plotting In P vs In pVo . The radius x = a;o ~ 0.01 co, 
corresponding to the minimal resolution of numerical simula- 
tions, marks the "inner" region (I). The "outer" region (II) 
denotes the halo up to its physical radius, the virial radius 
X = Co = 15. As a comparison we show a line with slope 1.18 
(thick grey line) that would correspond to the poytropic rela- 
tion with n ~ 5.5. Notice how the NFW halo approximately 
fits this relation, except near the center where the density 
profile becomes cuspy. 



VII. DISCUSSION AND CONCLUSSION 

The fact that NFW halos asymptoticaly comply with a 
polytropic relation with n « 5.5 is quite significant, since 
stellar polytropes characterized by (|45|l are the equilib- 
rium state associated with the entropy functional in the 
non-extensive entropy formalism derived by Tsallis and 
coworkers 31, 32, 33]. In its application to self- 
gravitating coUisionless systems this formalism is char- 
acterized by the free parameter q = {2n — l)/{2n — 3), 
so that the isothermal sphere (equilibrium state for the 
usual Boltzmann-Gibbs entropy functional) follows in 
the "extensive entropy" limit n ^ oo (or g — > 1). Assum- 
ing Tsallis theory to be correct, the empiric verification 
(see Figure^ that NFW halos outside the "inner" core 
satisfy a polytropic relation might indicate that in this 
"outer" region the NFW numerical simulations yield self- 
gravitating configurations that approach an equilibrium 
state characterized by the Tsallis parameter q « 1.25. 
However, while the central cusps in the density profile 
that are predicted by NFW simulations seem to be at 
odds with observations 0, [23| , there is no conflict be- 
tween these observations and the 1/x^ scaling of the 
NFW density profile outside the core region (as well as 



the rotation velocity profile from H4()(l '). Although the 
issue of the cuspy cores is still controversial, if galac- 
tic halos seem to exhibit fiat density cores, their profiles 
could be adjusted to stellar polytropes and this might 
be helpful in providing a better empirical verification of 
Tsallis' formalism. However, this idea must be handled 
with due case, since stellar polytropes follow from an 
isotropic velocity distribution, while galactic halos with 
such distributions could be unrealistic. 

As pointed out before, the density profile of NFW ha- 
los diverges at the center. Apparently this issue has not 
bothered astrophysicists, since (as mentioned before) the 
cuspy cores of NFW numerical simulations are meant to 
show a density scaling of 1/x near the center and these 
simulations cannot resolve distances to the halo center 
smaller than 1 % of the virial radius . One way to 
deal with this unphysical feature, leading to a better de- 
scription of these halos, would be to perform a smooth 
matching between NFW spacetimes and a small central 
region with a regular density profile. An adequate radius 
for this "inner" region could be the minimal resolution 
scale in numerical simulations {x — xq 0.01 cq). An- 
other improvement could be a smooth matching of the 
NFW spacetime to a Schwarzschild vacuum exterior at 
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the virial radius x = cq, which is the physical radius of 
the halo. One of the matching conditions in this lat- 
ter case would be 'P(co) = 0, implying a different choice 
of the integration constant C in H43|) . Another neces- 
sary improvement is the study of the anisotropic cases 
for which a ^ 0. 

We have constructed the spacetime corresponding to 
post-Newtonian generalizations suitable to NFW halos. 
Although we have presented only the idealized case with 
isotropic pressure, the methodology that we followed here 
can be applied, in principle, to any Newtonian model 
of galactic halos. We believe that it is necessary to 
study galactic halo models (NFW, as well as other em- 
piric or theoretical models) within a wider framework 
including the usual thermodynamics of self-gravitation 
systems ,2^ ,2^ , as well as alternative approaches such 



as Tsallis' formalism Such an improve- 

ment and extension of the present study of NFW halos 
are being pursued elsewhere [23 |. 
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